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Topological crystalline insulators (TCI) possess electronic states protected by crystal symmetries, 
rather than time-reversal symmetry. We show that the transition metal oxides with heavy transition 
metals are able to support nontrivial band topology resulting from mirror symmetry of the lattice. 
As an example, we consider pyrochlore oxides of the form A2M2O7. As a function of spin-orbit 
coupling strength, we find two Z2 topological insulator phases can be distinguished from each other 
by their mirror Chern numbers, each of which indicates a different TCI. We also derive an effective 
k ■ p Hamiltonian, similar to the model introduced for Pbi_ x Sn x Te, and discuss the effect of an on- 
site Hubbard interaction on the topological crystalline insulator phase using slave-rotor mean-field 
theory, which predicts new classes of topological quantum spin liquids. 
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Introduction.— A large class of materials has been pre- 
dicted to possesses a nontrivial Z 2 topological classifica- 
tion [1-5] resulting from an odd number of spin-orbit in- 
duced band inversions in the bulk of the material. These 
materials have gapless surface modes protected by time- 
reversal symmetry. However, crystal symmetries can also 
impose topological features, and therefore lead to an ad- 
ditional source of topological protection in materials pos- 
sessing these symmetries [6-9]. One class of such in- 
sulators is topological crystalline insulators (TCI) with 
surface states protected by either point [10] or mirror 
symmetry [11]. The recent discovery of TCI in the nar- 
row band semiconductor Pbi_ x Sn x Te/Se has paved the 
way to a new class of topological materials beyond the 
Zi topological insulators (TI) [12-14]. Thus, it is impor- 
tant to extend the domain of candidate TCI to include 
other materials, including those that may spontanously 
break time reversal symmetry (TRS) through magnetic 
ordering. 




FIG. 1. (Color online) Left: Pyrochlore lattice. Transition 
metal ions (large balls) are located on the corner shared tetra- 
hedra, each surrounded by an oxygen cage (red balls). Rear- 
earth elements fill spaces between cages. Two mirror planes 
are also shown. Center: Brillouin zone of face centered cu- 
bic lattice with mirror plane. Capital letters denote the high 
symmetry points. Left: a representation of unit cell with for 
sites and two possible magnetic orderings. On top ordering 
preserve mirror symmetry and on bottom ordering is of all- 
in-all-out which breaks mirror symmetry. 



In this work, we show that transition metal oxides 
(TMO) with hd orbitals are candidate materials for re- 
alizing a TCI, which extends the known list of potential 
TCI materials to those where interaction effects could 
play a larger role in driving qualitatively new phases and 
phenomena. Pyrochlore iridates, such as A 2 Ir 2 07, where 
A is a rear-earth element, have been shown to possess 
a nontrivial Zi invariant [15-17] so long as the strength 
of direct d-d hopping relative to indirect hopping via the 
oxygen orbitals does not fall within a certain window [18] 
and time-reversal symmetry is not broken [19]. 

Here, we explicitly show that surface states exist even 
if TRS is broken by the application of a magnetic field 
or from local magnetic moments induced by interactions. 
We argue that these surface sates are protected by mirror 
symmetry, and originate from a nonzero mirror Chern 
number [7, 11]. This nontrivial band topology leads to a 
new topological Mott phase in the regime of intermediate 
Hubbard interaction where the charge degrees of freedom 
on the surface become gapped, but the TCI structure of 
the gapless surface spin modes modes remain intact (a 
type of topological quantum spin liquid): the topological 
crystalline Mott insulator (TCMI). 

Model— The transition metal oxide A2R2O7 has a py- 
rochlore lattice structure (see Fig.l) composed of cor- 
ner shared tetrahedra of Ir +4 ions (<i-shell filling f=l/2). 
Each Ir +4 is surrounded by an octahedral cage of oxy- 
gen which splits the atomic d-levels into a lower t^g 
and higher e g manifold; spin-orbit coupling further acts 
within the t^g states to create eigenstates of total angular 
momentum. A minimal tight-binding model describing 
hopping of electrons between transition ions is as follows, 

i <ij> 

where first term contains both onsite energy and lo- 
cal spin-orbit coupling with strength A, and the second 
term describe nearest-neighbor hopping between transi- 
tion metal ions [15, 16, 18]. Here, 7 is a collective index 
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including both t2 g (yz, zx, xy) orbitals and spins, T ,ij, 
Td.ij are matrices for the oxygen mediated and direct 
hopping integrals, respectively. For simplicity, we con- 
sider only oxygen mediated hopping of electrons which 
captures the correct physical picture over a wide range 
of direct d hopping values [18]. Outside of this range, the 
system is gapless and a Weyl semi-metal can result from 
interaction effects [18, 20]. Physical pressure may provide 
a route to tune between these two interesting regimes. 

Because of the cubic structure of the lattice (face cen- 
tered cubic with four point basis), the lattice is sym- 
metric with respect to a variety of mirror planes. Two 
of them are shown in Fig.l. Crystal surfaces such as 
(010) are symmetric about the mirror planes, and we will 
show they support topologically protected gapless sur- 
face modes. The mirror operator transforming orbitals 
and spin about the mirror plane (101) is M = R (g> U, 
where R transforms the orbitals in local coordinates of 
each rotated octahedra and U rotates spin by 180° about 
an axis normal to the mirror plane, 

( T\ o o o \ 

R= ° ° ° r 2 u= i(-ll\ (2) 
o o r 3 o ^2 V 1 1 / 

V O T4 O O ) 

where is a reflection matrix transforming local ti g or- 
bitals on sites i—l.A (see Fig.l), and o is a 3x3 zero 
matrix. The rj are given by 
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The transformation of spin under the matrix U is used 
to determine which pattern of magnetic orderings or 
magnetic field preserve the mirror symmetry about the 
(101) plane. Note the reflection under this mirror plane 
also takes a unit cell located at lattice vector R^ = 
(xi,yi,Zi) to Ri = (Zi,yi,Xi). One can verify the k- 
space Hamiltonian transforms as MHd(k X) k y) fc z )M _1 = 
Hd(k Zl k y ,k x ). Thus, on the mirror plane in Brillouin 
zone (see Fig.l), where k x = k z , the Hamiltonian is in- 
variant, i.e. [Hd, M]=0, and, therefore, at each momen- 
tum on this plane the Bloch eigenstates of Hd can be 
labeled by mirror eigenvalues ±i since M 2 = — 1 [7, 11]. 

Mirror Chern number.— For 0<A<2.8 (we set NN hop- 
ping integral t=l) the non-interacting model is metallic; 
larger values of A open a gap in the spectrum and re- 
sult in a Z 2 TI with indicies (1;000) [15]. This gap 
closes for A c «3.25 by forming nodes at the four equiv- 
alent L points in Brillouin zone. These nodes will again 
be gapped for A>A C , but inverted relative to A<A C . How- 
ever, the Z2 index remains unchanged because an even 
number of band inversions occured. (See Table. I.) This 
gap closing phenomenon persists even in the presence of 
interactions [15]. In following, we will show these seem- 
ingly identical TI phases can be distinguished by a mirror 



TABLE I. The parity invariants at time reversal invariant 
momenta (TRIM) of bulk Bloch wavefunctions <5(r a ) at dif- 
ferent values of A. 
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Chern number. A similar gap closing at L points also oc- 
curs in Pbi_ x Sn x Te, separating a trivial insulator from 
a TCI [11]. 

The fact that all bands on the mirror plane can be la- 
beled by mirror eigenvalue ±« means that each band is 
effectively spin polarized. We have calculated the av- 
erage values of the spin operator <S> for each band 
and found for mirror related bands <S>oc±^(x — z), 
which is perpendicular to mirror plane. We also calcu- 
lated the Berry curvature f2(k) = V x A, where A = 
« ^ n (u„(k)|V|u„(k)) is the Berry connection (summed 
over all occupied bands) on the mirror plane in the Bril- 
louin zone (see central panel in Fig.l). In Fig.2a-b we 
plot fi (k) for occupied bands with mirror eigenvalues 
±i for A<A C and A>A C , respectively. The path in k- 
space is chosen along the mirror Brillouin zone bound- 
ary including L\ — (n, 7r,7r) and L2 = (tt, — 7r, 7r). The 
main contribution to fi(k) comes from k points around 
L. The plots clearly reveal that f2 +l (k) — Sl~ 4 (k) changes 
sign upon gap closing. For each polarization we calcu- 
lated the Chern number [7], um — { n +i ~ n—i)/2, and 
found that for A < A c : n + i=+l and n_j=-l yielding 
riAf=+l and for A > A c : ?im=-1. The nonzero value 
of mirror Chern number clearly shows that the mirror 
symmetry gives rise to nontrivial band topology which 
cannot be adiabatically connected to a trivial band insu- 
lator. Moreover, it shows that the Z2 TI phases around 
the gap closing point can be further labeled by UAf=±l. 

This behavior persists even in the presence of time 
reversal breaking perturbations. Time reversal symme- 
try can be violated by either applying a magnetic field 
or from local magnetic moments on the sites of the 
lattice. As an example, we add H' = Yli ■ Sj to 
Eq.(l). In order to understand which patterns of local 
moments or magnetic field preserve the mirror symme- 
try, we note that under reflection about mirror plane: 
o~ x — > —a z ,<j y — > —a y ,o- z — > — <j x . Hence, H' is invari- 
ant if Bi=±B(-l,0,l), or any magnetic ordering which 
is perpendicular to mirror plane. Two configurations of 
magnetic orderings are shown in Fig.l. The ordering on 
the top has mirror symmetry, while the all-in-all-out con- 
figuration breaks the mirror symmetry. In principle, such 
orderings can be stabilized by interactions [18, 20] 

First we consider a configuration of moments that pre- 
serve the mirror symmetry. The mirror eigenvalues are 
still well defined and can be used to label the Bloch states. 



3 



a 





i 


(b) / 


\ / 




i-i- 



-2 




/ +i \ 






Li 




U 




Li L 2 


(c) 








0.3 


(d) 










-i 


A 















+i 




-0.3 










Li 




U 




L, L 2 



FIG. 2. (Color online) Berry curvature fi(k) along the bound- 
ary of the mirror plane in Brillouin zone. L\ — (n, ir, n) and 
1/2 = (tt, — 7r, 7r) are high symmetry points of 3D Brillouin 
zone residing on the boundary of the mirror plane in k-space. 
In (a)-(b) the time reversal symmetry is preserved and A<A C 
in (a) and A>A C in (b). In (c)-(d) the time reversal symme- 
try is broken by the magnetic ordering shown in Fig.l which 
preserves and breaks the mirror symmetry in (c) and (d), re- 
spectively. 



In Fig. 2c we plot the corresponding Berry curvatures. Al- 
though these values of Berry curvature look asymmetric 
between states with mirror eigenvalues +i and — i , the 
Berry curvatures integrated over entire mirror Brillouin 
zone are the same up to a minus sign, namely n + i = —n—i 
yielding a nonzero integer value nA/=-l. This identifi- 
cation shows that the band topology is constrained by 
mirror symmetry not by TRS. If the mirror symmetry is 
broken, say by an all-in-all-out ordering, the Bloch sates 
are no longer eigenstates of the mirror operator. The cor- 
responding Berry curvature is shown in Fig. 2d. Due to 
its antisymmetric nature around each L point, the Berry 
phase of each L point and thus the total Chern number is 
zero. While the two magnetic orderings give rise to zero 
total Chern number, the mirror symmetric one yields an 
integer value of mirror Chern number. 

Surface states.— The bulk-boundary correspondence 
implies the surface between two insulators with differ- 
ent bulk band topology carries surface modes [3, 4]. In 
order to confirm that the surface of the pyrochlore lat- 
tice supports gapless surface states we consider a slab 
geometry of the lattice along the [010] direction. This 
slab, as shown in Fig.l, is symmetric about two mirror 
planes: One is (101) and another one is (101) with an off- 
set of 1/4 (or 3/4) of lattice spacing. We diagonalize the 
Hamiltonian in this geometry and plot the surface states 
along the high symmetry points of the projected Bril- 
louin zone. The results for different sets of parameters 




FIG. 3. (Color online) The band structure of the tight-binding 
model in Eq.(l) along the high symmetry points of projected 
Brillouin zone as shown in inset of (a). In panels (a)-(b) the 
time reversal symmetry is preserved and A=3 (a) A=4 (b). In 
panels (c)-(e) the time reversal symmetry is broken, but the 
mirror symmetry is preserved: (c) magnetic field B=B(1,0,1), 
(d) magnetic field B=B(-1,0,1) and (e) magnetic ordering. In 
panel (f ) both time reversal and mirror symmetries are broken 
by a magnetic ordering of all-in-all-out type. 



are shown in Fig. 3. For a time reversal symmetric system 
there are Dirac nodes right at the projected time rever- 
sal invariant momenta (TRIM) M and X (see Fig.3a-b). 
By increasing spin-orbit coupling across the gap closing 
point at A c the overall features of the surface states re- 
main unchanged except that the area of Fermi surface 
enclosed by Dirac nodes at M change and its electron or 
hole like character on the two surfaces of the slab change. 

Now we break the time reversal symmetry by adding 
H' . We first consider the cases in which the perturba- 
tion respects the mirror symmetry. The results for uni- 
form magnetic field B!=B(1,0,1) and B 2 =B(-1,0,1) are 
shown in Fig.3(c)-(d), respectively. Note that Bi breaks 
the symmetry about mirror plane (101) but preserves the 
mirror symmetry about plane (101), while field B 2 does 
opposite. The Kramers degeneracy is lifted and Dirac 
nodes become gapped at M and X point. However, ad- 
jacent to X there are still some crossing surface modes. 
Mirror planes in Brillouin zone, such as those shown in 
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Fig.l, define mirror lines in the projected Brillouin zone. 
For example, for k points on the line YX, where k x =k z , 
the eigenstates can be labeled by mirror eigenvalues. The 
same arguments hold for line MX. Thus, surface states 
can also be labeled by their corresponding mirror eigen- 
values. The crossing of the surface states along the MX 
line is protected by mirror symmetry about (101) plane 
(see Fig. 3c), and the crossing along YX line is protected 
by mirror symmetry about (101) plane (see Fig. 3d). Note 
that away from these lines the degeneracy is lifted. We 
also calculated the surface states in the presence of a 
pattern of magnetic ordering which preserves symmetry 
about the mirror plane (101). The result is shown in 
Fig.3e. The surface modes cross each other along YX 
line as expected. This identification clearly shows that 
the crossing surface modes are associated with mirror 
symmetry, independent of the details of a time reversal 
breaking perturbation. 

Finally, we break symmetry with respect to both mir- 
ror planes by considering an all-in-all-out magnetic con- 
figuration. The band structure is shown in Fig.3f. It is 
seen that the degeneracy is lifted and there is no crossing 
along either MX or YX. We also explored many other 
cases (not shown here). All mirror symmetry breaking 
perturbations remove the band crossing in the gap. This 
finding further verifies that the surface crossing discussed 
in the previous paragraphs are associated with mirror 
symmetry. 

kp Hamiltonian.— It is instructive to give a k • p ex- 
pansion of Eq.(l) near the L point. We first rotate 
{k x ,k y ,k z ) into a new orthogonal basis (^1,^2,^3) in 
which k\ is along the YL line and k% is perpendicular 
to mirror plane. Thus, the mirror plane corresponds to 
&3 = 0. Right at the L point Eq.(l) in k-space becomes 
block diagonal with Bloch eigenstates localized on either 
site 1 (or sites 234) [see unit cell in Fig.l]. This effectively 
subdivides the lattice sites in two sets and defines a basis 
to expand Eq.(l) around the L point. We label them 
as tpA and ipB> respectively. With this identification the 
effective Hamiltonian near the L point is, 

H ±i = mr* ± vi • kr* +v 2 -kr !/ , (4) 

where t z — ±1 corresponds to tpA an d ipB , respec- 
tively, and k = (fci,fc 2 ). Eq.(4), and therefore its un- 
derlying physics, is not dissimilar from the model intro- 
duced to describe the insulator phases of semiconductor 
Pbi_ x Sn x Te upon doping [11]. In fact, in both mod- 
els the character of conduction/ valence bands- (i-orbital 
V'a/Vb in °ur model and p-orbital cation/onion in the 
Pbi_ x Sn x Te model-gets switched at L point by tuning 
spin-orbit coupling A and Pb doping, respectively. The 
band inversion changes the sign of the mass term m in 
Eq.(4), which in turn changes the Chern number for ±i 
mirror eigenstates by ±1. Since the gap closing occurs at 
two equivalent L points (related by 180° rotation about 
[101] axis) on the mirror plane, the total change of Chern 



number for each ±i states will be ±2. This change of 
Chern number is consistent with Berry curvatures calcu- 
lated in Fig. 2. There is, however, a significant difference 
between our model and the model for Pbi_ x Sn x Te. In 
latter model, the gap closing at L points signals a topo- 
logical phase transition between a trivial insulator PbTe 
and topological insulator SnTe with mirror Chern num- 
ber tim = 2. But in our d-orbital model the topological 
phase transition occurs between two topological phases 
with different mirror Chern number: tim = — 1 before 
and riM = +1 after a gap closing. This distinction be- 
tween topological insulator phases has some implications 
for the surface states which can potentially be explored 
in experiment. On the mirror line XYX there is exactly 
one crossing surface mode due to n^i = +1, while for 
Pbi_ x Sn x Te there are two crossings. 

Topological crystalline Mott insulator.— In d-orbital 
models the electron correlations are also strong and 
might even lead to new topological phases, such as the 
topological Mott insulator (TMI) [15]. The TMI phase 
can be obtained through a spin-charge separation inher- 
ently included in the slave-rotor decomposition [21, 22]: 
Cj y = e l6 i fjj. The mean-field theory can be written as 
H = Hf + Hg, where Hf and Hg describe, respectively, 
spinon and rotor parts of the model, and are related to 
each other via mean- field parameters Qf — (e~ l ( 0i ~ ^) 

and Q g = T ^ ft,afji'a') which are determined 

self-consistently as a function of Hubbard interaction 
[15, 16]. The Hamiltonian of the electrons is described 
by Eq.(l) with e^ 7 — > /j 7 and 7y — > QfTij. Thus, the 
electron correlations systematically renormalize the band 
width of the system by decreasing Qf. Mott physics oc- 
curs when quasiparticle weight Z =(e lS ) vanishes for suf- 
ficiently strong Hubbard interactions, where the bosonic 
rotor excitations are gapped and become uncondensed. 
Beyond this limit the low energy excitations are described 
solely by a spinon Hamiltonian, Hf, which may have 
topological band structure. We calculated the mirror 
Chern number in the boson uncondensed phase and es- 
tablished that um = — 1- This already means that the 
previous topological crystalline insulator phase at weak 
interacting limit will transit into a Mott analogue-the 
topological crystalline Mott insulator (TCMI)-which is 
similar in physics to a topological Mott insulator phase 
inherited from the nontrivial band topology of a nonin- 
teracting TI only with gapless spin modes protected by 
the mirror symmetry rather than TRS. 

Summary. —In this work we discussed the possible real- 
ization of topological crystalline insulators in transition 
metal oxides with 5d orbitals and mirror symmetry, and 
shown electron interaction effects can drive a novel TCMI 
phase. We thank Liang Fu for helpful discussions and 
we gratefully acknowledge financial support from ARO 
Grant No. W911NF-09-1-0527, W911NF-12-1-0573, and 
NSF Grant No. DMR-0955778. 
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